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ABSTRACT 


The  existence  of  two  motive  forces  on  a  Crookes  radiometer  has  complicated  the 
investigation  of  either  force  independently.  The  thermal  creep  shear  force,  in  particular, 
has  been  subject  to  differing  interpretations  of  the  direction  in  which  it  acts  and  its  order 
of  magnitude.  A  horizontal  vane  radiometer  design  is  provided,  which  isolates  the 
thermal  creep  shear  force.  The  horizontal  vane  radiometer  is  explored  through 
experiment,  kinetic  theory,  and  the  Direct  Simulation  Monte  Carlo  (DSMC)  method.  The 
qualitative  agreement  between  the  three  methods  of  investigation  is  good.  The 
quantitative  agreement  between  the  three  methods  of  investigation  is  better  than  an  order 
of  magnitude  in  the  cases  examined.  The  thermal  creep  force  is  found  to  act  from  the  hot 
side  to  the  cold  side  of  the  vane.  The  peak  in  the  radiometer’s  angular  speed  as  a  function 
of  pressure  is  found  to  be  explained  as  much  by  the  behavior  of  the  drag  force  as  by  the 
behavior  of  the  thermal  creep  force.  In  addition,  this  dissertation  provides  scaling  laws 
between  millimeter-scale  and  micron-scale  horizontal  vane  radiometers,  a  design  of  a 
microelectromechanical  system  (MEMS)  horizontal  vane  radiometer,  and  conceptual 
designs  of  two  MEMS  energy  harvesting  devices  that  exploit  the  thermal  creep  force. 
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I.  INTRODUCTION 


The  purpose  of  this  dissertation  is  twofold:  first,  to  confirm  and  improve  upon  the 
existing  understanding  of  the  thermal  creep  force,  and  second,  to  design 
Microelectromechanical  System  (MEMS)  energy  harvesting  devices  that  utilize  the 
thermal  creep  force.  Recent  research  using  numerical  methods  that  shows  the  force  acting 
in  a  direction  opposite  to  that  originally  found  by  Osbourne  Reynolds  motivates  the  need 
to  confirm  the  existing  understanding  of  the  thermal  creep  force  [1],  [2].  The  fact  the 
solution  to  the  Boltzmann  equation  with  the  Chapman-Enskog  correction  for  thermal 
creep  force  has  been  reported  to  significantly  over-predict  experimental  results  and  has 
not  been  fully  accepted  by  the  community  motivates  the  need  to  improve  on  the  existing 
understanding  of  the  thermal  creep  force  [3]. 

This  chapter  provides  a  historic  overview  of  the  Crookes  radiometer,  which  led  to 
the  discovery  of  the  thermal  creep  force.  Chapter  II  presents  the  results  of  an  experiment 
on  a  novel  horizontal  vane  radiometer  designed  to  isolate  the  thermal  creep  force  from 
the  other  motive  force  on  the  Crookes  radiometer.  These  results  confirm  Reynolds’s 
report  of  the  direction  in  which  the  thermal  creep  force  acts.  Chapter  III  considers  several 
edge  and  kinetic  theory  corrections  to  the  thermal  creep  force  result  given  by  Scandurra 
et  al.  [4].  These  corrections  allow  the  result,  when  combined  with  the  drag  force,  to 
provide  a  reasonable  prediction  of  the  angular  speeds  reported  in  Chapter  II.  In  Chapter 
IV,  the  results  of  numerical  simulations  of  the  experiment  in  Chapter  II  are  presented. 
The  purpose  of  the  simulations  is  to  provide  direct  values  of  the  thermal  creep  and  drag 
forces.  These  results  lend  credence  to  the  adjustments  made  in  Chapter  III.  In  Chapter  V, 
a  comparison  of  experimental,  theoretical,  and  simulation  results  is  given.  Finally, 
Chapter  VI  provides  a  scaling  law  for  the  thermal  creep  force  on  MEMS  vanes  and  two 
designs  of  energy  harvesting  devices. 

A.  THERMAL  CREEP 

Thermal  transpiration,  or  thermal  creep,  is  the  flow  of  a  gas  over  a  surface  due  to 
a  temperature  gradient  in  the  gas  parallel  to  the  surface.  The  force  on  the  gas  is  equal  and 
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opposite  to  the  force  on  the  surface.  The  force  exerted  on  a  surface  by  a  gas  with  a 
temperature  gradient  is  thus  called  the  thermal  creep  force.  The  force  is  most  famously 
demonstrated  by  the  Crookes  radiometer  (see  Figure  1). 


Figure  1.  A  Crookes  Radiometer 


The  figure  on  the  right  is  an  edge  view  of  the  vane  and  not  to  scale.  In  commercial  Crookes  radiometers, 
typical  vanes  are  1.5  cm  x  1.5  cm  and  thickness  is  usually  on  the  order  of  0.1  mm.  The  thermal  creep  force 
acts  on  the  edge  parallel  to  the  temperature  gradient.  The  Einstein  effect  acts  on  the  surfaces  perpendicular 
to  the  temperature  gradient. 

William  Crookes  published  the  atomic  weight  of  Thallium  in  1873.  He  weighed 
Thallium  salts  at  atmospheric  pressure  and  reduced  pressure  to  quantify  the  buoyant  force 
on  his  sample  and  determine  what  the  weight  would  be  in  a  perfect  vacuum.  He  found 
that  the  weight  measured  at  reduced  pressures  was  sensitive  to  the  light  source  he  used  to 
illuminate  the  scale.  At  the  time  it  was  believed  that  a  hot  object  would  weigh  less  than  a 
cold  object  of  similar  mass  because  of  upward  air  flow  near  the  hot  object  or 
condensation  on  the  cold  object  [5].  Crookes  recognized  these  effects  should  decrease 
with  decreasing  pressure,  but  the  effect  he  was  observing  did  not.  He  designed  the 
Crookes  radiometer  and  other  experiments  to  study  the  effect.  In  1874  he  ascribed  the 
observed  motion  to  repulsion  by  radiant  heat  which  he  theorized  is  present  at  all  pressures 
but  masked  by  air  currents  at  higher  pressures  [6].  This  is  not  the  correct  explanation  of 
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the  motion.  Though  radiation  pressure  does  exist,  it  produces  a  small  force  and  would 
tend  to  move  the  vanes  with  the  black  side  leading  contrary  to  the  observed  motion. 

The  same  year  Crookes  published  his  results  and  analysis,  Osbourne  Reynolds 
published  an  analysis  ascribing  the  effects  to  the  transfer  of  momentum  between  the 
surface  and  gas  molecules  during  evaporation,  condensation,  and  heat  transfer  [7]. 
Osbourne  Reynolds  further  explains  that  an  evaporating  molecule  leaves  the  surface  with 
a  certain  momentum  and  the  surface  gains  an  equal  and  opposite  momentum.  A  molecule 
undergoing  condensation  arrives  with  a  certain  momentum  which  is  transferred  to  the 
surface.  A  molecule  that  transfers  heat  with  a  surface  during  a  collision  departs  with  a 
different  momentum  than  with  which  it  arrived  the  balance  of  which  is  conserved  by 
transfer  to  or  from  the  surface. 

Reynolds  continued  to  study  the  transfer  of  momentum  between  a  surface  and  gas 
during  heat  transfer  by  experimenting  with  a  porous  plate  with  its  two  sides  held  at 
different  temperatures  [8].  He  was  able  to  conduct  his  experiments  at  atmospheric 
pressure  because  the  ratio  of  a  characteristic  length  of  the  pores  to  the  mean  free  path  of 
air  at  atmospheric  pressure  is  similar  to  the  ratio  of  the  characteristic  length  of  a  vane  in  a 
Crookes  radiometer  to  the  mean  free  path  of  air  at  the  pressures  where  motion  can  be 
observed.  He  found  that  gas  would  flow  through  the  pores  from  the  cold  side  of  the  plate 
to  the  hot  side  of  the  plate  and  named  this  flow  “thermal  transpiration”  in  1879. 
Molecules  entering  the  pores  from  the  cold  side  of  the  plate  tend  to  gain  heat  and 
therefore  velocity  as  they  travel  through  the  pore  while  molecules  entering  the  pores  from 
the  hot  side  of  the  plate  tend  to  lose  heat  and  therefore  velocity  as  they  travel  through  the 
pore.  These  effects  combine  for  a  net  velocity  and  flow  from  the  cold  side  to  the  hot  side. 
This  same  flow  occurs  over  the  edge  of  the  vanes  in  a  Crookes  radiometer  creating  a 
thermal  creep  force  on  the  vane  directed  from  the  hot  side  to  the  cold  side. 

B.  THE  ENERGY  ACCOMMODATION  COEFFICIENT 

Also  in  1879,  but  after  reading  Reynolds’  paper,  James  Maxwell  further  refined 
the  theory  of  thermal  creep  by  considering  the  details  of  the  collisions  between  gas 
molecules  and  the  surface  [9].  Maxwell  considered  two  types  of  collisions.  Reflected 
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gases  behave  as  if  the  surface  was  perfectly  smooth  and  elastic.  No  energy  is  transferred 
between  a  reflected  gas  and  the  surface,  and  the  tangential  components  of  the  molecules’ 
velocities  are  unchanged  by  the  collision.  Reflected  gas  does  not  contribute  to  thermal 
creep.  Absorbed-and-evaporated  gases  behave  as  if  they  were  “entangled  in  the  stratum 
of  spheres”  that  Maxwell  imagined  coving  the  surface.  Energy  is  transferred  between  an 
absorbed-and-evaporated  gas  and  the  surface.  Absorbed-and-evaporated  gas  molecules 
leave  the  surface  in  a  random  direction  with  a  speed  equivalent  to  the  surface 
temperature.  Absorbed-and-evaporated  gas  is  responsible  for  thermal  creep.  The 
magnitude  of  the  thermal  creep  effect  is  thus  related  to  the  ratio  of  the  number  of  gas 
molecules  that  are  absorbed-and-evaporated  to  the  total  number  of  gas-surface  collisions 
in  a  sample. 


Today,  the  extent  to  which  energy  is  transferred  in  a  collision  between  a  gas 
molecule  and  a  surface  is  expressed  as  the  accommodation  coefficient  a.  An 
accommodation  coefficient  of  zero  is  equivalent  to  a  perfectly  reflected  gas.  An 
accommodation  coefficient  of  one  is  equivalent  to  a  perfectly  absorbed-and-evaporated 
gas.  Other  models  use  separate  accommodation  coefficients  for  energy,  normal 
momentum,  and  tangential  momentum  calculations;  but  the  Maxwell  model  implies  those 
values  are  the  same.  The  accommodation  coefficients  of  real  gases  and  surfaces  are  a 
function  of  the  colliding  molecules  and  the  temperature.  The  modified  Boule  formula 
estimates  the  accommodation  for  clean  surfaces  below  500  K,  where  p  is  the  ratio  of  the 
molecular  weight  of  the  gas  to  the  molecular  weight  of  the  surface. 


a  = 


2.4// 

(i AT 


(u) 


Surfaces  absorb  layers  of  molecules  when  exposed  to  air  [10].  A  gas  molecule  colliding 
with  a  “dirty”  surface  will  thus  collide  with  a  molecule  of  a  similar  molecular  weight  (p 
equals  one)  regardless  of  the  material  composition  of  the  underlying  surface.  From 
equation  (1.1),  the  accommodation  coefficient  equals  its  maximum  value  of  0.6  when  p 
equals  one.  Thus,  0.6  is  a  reasonable  estimate  for  the  energy  accommodation  coefficient 
of  a  “dirty”  surface  in  air.  Nevertheless,  it  should  be  noted  accommodation  values  greater 
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than  the  maximum  value  of  the  modified  Boule  formula  have  been  observed  and  are 
commonly  assumed  in  other  studies. 

C.  THE  EINSTEIN  EFFECT 

Thermal  creep  is  a  shear  force  acting  on  the  edges  of  the  vanes  in  a  Crookes 
radiometer.  In  1924,  Albert  Einstein  identified  a  normal  force  acting  on  the  perimeter  of 
the  vertical  faces  as  well  [11].  He  imagined  a  chamber  of  rarified  gas  with  a  temperature 
gradient  and  no  net  mass  flow.  A  flat  surface  smaller  than  a  mean  free  path  of  the  gas 
placed  in  the  chamber  normal  to  the  temperature  gradient  would  experience  a  force 
toward  the  cold  side  as  the  same  number  of  molecules  per  unit  time  would  collide  with  it 
from  the  cold  side  as  from  the  hot  side  because  of  the  no  net  mass  flow  requirement  while 
the  particles  coming  from  the  hot  side  would  have  a  higher  velocity.  A  flat  surface  larger 
than  a  mean  free  path  of  the  gas  in  the  same  chamber  oriented  in  the  same  way  would  not 
have  a  net  force  on  its  center  once  the  pressure  had  time  to  equalize  and  reestablished  the 
no  mass  flow  requirement.  There  would  nevertheless  be  a  transition  region  between  the 
center  of  the  surface  where  no  net  force  is  experienced  and  the  region  beyond  the  surface 
where  a  small  surface  placed  in  the  chamber  would  experience  a  force.  This  transition 
occurs  over  a  mean  free  path  wide  perimeter  of  the  surface.  The  vanes  of  the  Crookes 
radiometer  have  a  net  normal  force  in  this  perimeter  directed  from  the  hot  side  to  the  cold 
side.  A  mathematical  development  of  the  Einstein  effect  is  included  in  Appendix  A. 

Thermal  creep  and  the  Einstein  effect  provide  the  theoretical  bases  for  the  two 
motive  forces  on  the  Crookes  radiometer.  Subsequent  work  on  the  topic  has  focused  on 
mathematical  refinement  of  the  theories  and  novel  radiometer  designs  [3].  One  alternative 
design  is  the  Hettner  radiometer  (see  Figure  2).  In  1924,  M.  Czerny  and  G.  Hettner 
designed  an  experiment  to  measure  the  shear  pressure  or  force  on  a  plate  with  a  uniform 
temperature  from  the  thermal  creep  caused  by  a  parallel  plate  with  a  temperature  gradient 
[12].  This  design  does  not  experience  the  Einstein  effect  because  the  plate  face  is  parallel 
to  the  temperature  gradient. 
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Figure  2.  Hettner  Radiometer 
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Plate  1  has  a  temperature  gradient  and  plate  2  has  a  uniform  temperature.  The 
temperature  gradient  establishes  thermal  creep  (arrows  with  fletchings).  Plate  1 
experiences  thermal  creep  force  while  plate  2  experiences  drag  from  the  thermal  creep. 
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II.  A  HORIZONTAL  VANE  RADIOMETER  EXPERIMENT 


We  designed  a  horizontal  vane  radiometer  (see  Figure  3)  that,  like  the  Hettner 
radiometer,  has  the  temperature  gradient  parallel  to  the  vane  face  but,  like  the  Crookes 
radiometer,  is  free  to  rotate  around  a  spindle  [13]. 


Figure  3.  The  Crookes  and  Horizontal  Vane  Radiometers 


(a)  The  Crookes  radiometer  experiences  thermal  creep  forces  and  the  Einstein  effect. 

(b)  The  horizontal  vane  radiometer  experiences  thermal  creep  forces  but  no  Einstein 
effect. 


The  angular  speeds  of  horizontal  vane  radiometers  were  observed  over  a  range  of 
pressures  with  mean  free  paths  from  5  mm  to  0.2  mm.  Two  types  of  horizontal  vane 
radiometers  were  constructed.  The  narrow  vane  radiometer  had  8  mm  by  16  mm  vanes. 
The  wide  vane  radiometer  had  16  mm  by  16  mm  vanes.  The  vanes  were  constructed  from 
high  gloss  photo  paper.  Half  of  the  non-glossy  side  was  printed  with  black  LaserJet  ink. 
Two  pieces  of  paper  cut  into  the  shape  of  the  vanes  were  glued  together  glossy  side  to 
glossy  side  so  that  both  the  upper  and  lower  horizontal  surfaces  would  have  a  white  side 
and  a  black  side.  The  vanes  on  both  radiometers  were  connected  to  the  spindle  with  4  mm 
stems. 

A  high  intensity  radiant  heat  projector  was  used  to  illuminate  the  radiometers.  The 
temperature  profile  that  developed  on  the  vanes  when  illuminated  was  observed  with  a 
FLIR  SC8200  camera  outside  of  the  vacuum  chamber  and  a  FLIR  ThermoVision  A20 


7 


camera  inside  of  the  vacuum  chamber.  The  temperature  profile  was  not  observed  to  be 
sensitive  to  pressure.  Under  illumination,  the  vanes  developed  a  9  K  temperature 
difference  between  the  black  and  the  white  side  with  the  transition  in  temperature  mostly 
occurring  in  a  sharp  gradient  in  the  center  of  the  vane  with  a  characteristic  length  of  3.5 
mm.  The  radiometers  were  tested  in  a  18-inch  diameter  12-inch-high  Pyrex  bell  jar  on  a 
highly  reflective  test  stand  so  that  both  the  upper  and  lower  horizontal  surfaces  would  be 
illuminated.  Figure  4  reports  the  results  of  these  experiments  at  the  pressures  tested.  [14] 
shows  video  of  horizontal  vane  radiometers  rotating  at  10  mTorr. 


Figure  4.  Angular  Speed  (Experiment) 
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The  narrow  vanes  are  8  mm  by  16  mm  and  the  wide  vanes  are  16  mm  by  16  mm.  Both 
radiometers  have  four  vanes.  Vanes  of  both  radiometers  have  a  4  mm  stem  connecting 
them  to  the  spindle  axis.  The  temperature  difference  on  the  vanes  is  9  K.  The 
characteristic  length  of  the  temperature  gradient  is  3.5  mm.  The  observed  motion  is  with 
the  white  (cold)  side  leading.  Angular  speed  is  reported  in  revolutions  per  minute 
(RPMs). 
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III.  THE  HORIZONTAL  VANE  RADIOMETER  THEORY 


The  theory  of  thermal  creep  explains  why  the  horizontal  vane  radiometer  rotates, 
but  an  understanding  of  why  the  angular  speed  peaks  and  is  a  function  of  the  vanes’ 
widths  requires  careful  examination  of  the  details  of  both  the  drag  force  and  thermal 
creep  shear  force  on  the  vanes.  We  follow  the  development  of  Wolfe  et  al.  [15]. 


A.  THERMAL  CREEP 


For  an  angular  speed  of  10  radians  per  second,  the  furthest  edges  of  the  vanes  will 
move  a  distance  between  2  pm  and  80  nm  during  a  mean  free  time  between  collisions  for 
the  pressures  in  our  experiment.  Because  these  distances  are  small  relative  to  the  widths 
of  the  vanes  in  our  experiment,  a  temperature  gradient  in  the  gas  over  the  vane  created  by 
collisions  between  the  gas  and  the  vane  will  remain  over  the  vane  despite  the  vane’s 
rotation  justifying  a  quasi-static  assumption  for  the  calculation  of  the  thermal  creep  force 
on  a  moving  radiometer  with  the  observed  velocities. 


Molecules  in  a  gas  with  a  temperature  gradient  have  a  Chapman-Enskog 
distribution  of  velocities  [16].  The  probability  density  function  /  is  the  sum  of  the 
Maxwell-Boltzmann  distribution  function  fo  and  the  Chapman-Enskog  correction/}. 

f  =  fo  +  fi  (3-1) 


The  Chapman-Enskog  correction/}  for  a  gas  with  a  temperature  gradient  is 


fi  =  Avy  --/?(v,2  +vj  +vz2)  /c 


fo  = 


f  e-p{vw,™i) 


A  = 


fi  = 

15 


m 


2  kBT 


\n m  dT 


32 nacsT  \  kBT  dy 


(3.2) 

(3.3) 

(3.4) 

(3.5) 
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where  v  is  the  velocity,  m  is  the  mass,  n  is  the  number  density,  <jcs  =  nd 2  is  the  hard- 

sphere  collision  cross  section  of  the  molecules  with  diameter  d,  kB  is  Boltzmann’s 
constant,  T  is  the  gas  temperature,  and  dT/dy  is  the  temperature  gradient  in  the  gas  [16]. 


The  net  thermal  creep  pressure  Aprc  on  a  surface  in  the  y-z  plane  with  the  gas  in 
the  positive  x  direction  is  determined  by  integrating  the  longitudinal  momentum  transfer 
of  the  gas  along  the  surface  of  the  plate  with  accommodation  coefficient  a  [4] . 

0  oo  oo 

A pTC  =  nm  J  dvx  J  dvz  J  dvyvyvxf  (3.6) 


The  Maxwellian  term  integrates  to  zero  but  the  Chapman-Enskog  term  survives 
providing  the  following  result. 


A Ptc 


15 


kB  dT 


64>/2  crcs  dy 


a 


(3.7) 


The  thermal  creep  force  FTC  is  this  pressure  multiplied  by  the  area  over  which  it 
acts  Arc-  If  LTc  is  the  length  of  a  rectangle  perpendicular  to  the  temperature  gradient  and 
Wtc  is  the  width  along  the  temperature  gradient,  we  have 


k'rc  AprcA-rc  —  Aprc  LrcW./c . 


(3.8) 


Assuming  the  temperature  gradient  in  the  gas  is  roughly  constant  in  the  area 
where  the  thermal  creep  force  acts,  the  temperature  gradient  term  can  be  replaced  by 


dT  AT 

^  WGgas  ’ 


(3.9) 


where  AT  is  the  temperature  difference  between  the  hot  and  the  cold  sides  of  the  surface 
and  Wcgas  is  the  width  of  the  temperature  gradient  in  the  gas. 


The  width  of  the  temperature  gradient  established  in  a  low-pressure  gas  by  a 
surface  with  a  temperature  gradient  is  wider  than  the  temperature  gradient  on  the  surface 
Wcsurf  by  a  length  equal  to  the  slip  length  [17],  or 

22  (3.10) 
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where  X  is  the  mean  free  path  of  the  gas  and  a  is  the  accommodation  coefficient.  The 
mean  free  path  is  related  to  the  particle  number  density  n  and  collision  cross  section  as 


The  physical  meaning  of  equation  (3.10)  can  be  understood  by  considering  the 
two  extreme  values  of  a.  When  the  accommodation  equals  zero,  the  width  of  the  gradient 
in  the  gas  is  infinite.  There  is  thus  no  gradient  when  there  is  no  energy  transfer  between 
the  surface  and  the  gas.  When  the  accommodation  equals  one,  the  width  of  the  gradient  in 
the  gas  is  the  width  of  the  gradient  on  the  surface  plus  two  mean  free  paths.  The  gradient 
thus  extends  one  mean  free  path  further  on  each  side  in  the  gas  than  in  the  surface. 

In  the  horizontal  vane  radiometer  experiment,  there  are  four  temperature  gradients 
in  the  plane  of  each  vane  (see  Figure  5). 


Figure  5.  Temperature  Gradients  on  a  Horizontal  Vane 
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The  white  (cold)  side  of  the  vane  is  assumed  to  have  a  temperature  similar  to  ambient 
temperature  and  the  grey  (hot)  side  is  warmer  than  ambient  temperature. 
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Temperature  gradient  #1  causes  the  force  FTo  which  causes  the  rotation  of  the 
radiometer  and  is  a  convenient  starting  point  for  the  calculation  of  the  total  thermal  creep 
force  on  the  vane.  We  will  treat  the  effects  of  the  other  temperature  gradients  as  apparent 
modifications  of  the  force  from  temperature  gradient  #1. 

Temperature  gradient  #2  creates  a  force  which  opposes  the  rotation.  The 
magnitude  of  this  force  is  less  than  the  magnitude  of  the  force  created  by  temperature 
gradient  #1  because  the  portion  of  temperature  gradient  #2  which  is  not  over  the  vane 
does  not  create  a  force,  i.e. 

FTC=xFTCl,  (3.11) 

where  x  is  a  numerical  factor  between  zero  and  one  representing  the  fraction  of  the  force 
from  temperature  gradient  #1  which  is  offset  by  the  force  from  temperature  gradient  #2. 


The  net  force  created  by  temperature  gradients  #3  and  #4  is  zero,  but  these 
gradients  decrease  the  magnitude  of  temperature  gradient  #1  near  the  edges.  This  has  the 
effect  of  decreasing  the  apparent  length  over  which  the  net  thermal  creep  pressure  from 
temperature  gradient  #1  acts.  The  length  correction  is  the  slip  length  times  the  numerical 
factor  [)  or 


L.J  =  LIcl-pX— (3.12) 
a 

The  width  over  which  temperature  gradient  #1  acts  is  the  smaller  of  either  the 
width  of  the  vane  Wy  or  the  width  of  the  temperature  gradient  in  the  gas.  Putting  this 
together  with  equations  (3.7)  to  (3.12)  gives  the  following  expression  for  the  net  thermal 
creep  force  on  the  two  sides  of  a  single  vane  in  the  horizontal  vane  radiometer 
experiment: 


k 

iv,  = — — —AT  aL'  min 


TC 


32V2 


(7, 


CS 


Wy  ' 

V  ^Ggas  J 


where  L  is  the  length  of  the  vane  with  the  slip  length  correction. 


(3.13) 


Figure  6  shows  the  net  thermal  creep  force  predicted  by  this  equation  for  the 
vanes  in  the  horizontal  vane  radiometer  experiment  assuming  numerical  factors  of 
v  =  0.5  and  p  =  0.5  for  equations  (3.11)  and  (3.12).  These  values  follow  from  the 
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assumption  that  the  gas  temperature  at  the  edge  is  equal  to  the  average  of  the  temperature 
of  the  hot  side  of  the  vane  and  the  ambient  temperature.  We  assumed  an  accommodation 
coefficient  of  a  =  0.6  .  An  even  higher  value  of  the  accommodation  would  result  in  the 
force  rising  more  quickly  as  a  function  of  pressure  and  asymptotically  approaching  a 
higher  value.  The  pressure  dependence  of  the  force  is  a  result  of  the  slip  length 
corrections  in  the  width  of  the  gradient  in  the  gas  and  the  apparent  length  of  the  vane. 


Figure  6.  Thermal  Creep  Force  (Theory) 


The  narrow  vanes  are  8  mm  by  16  mm  and  the  wide  vanes  are  16  mm  by  16  mm.  The 
temperature  difference  on  both  vanes  is  9  K.  The  characteristic  length  of  the  temperature 
gradient  on  the  surface  is  3.5  mm.  The  accommodation  coefficient  is  0.6.  At  60  mTorr  and 
greater  pressure  the  forces  are  similar  on  the  narrow  and  wide  vanes.  A  positive  force  is 
directed  from  the  hot  side  to  the  cold  side. 


B.  DRAG  FORCE 

The  dominant  drag  force  on  a  thin  plate  oriented  parallel  to  the  flow  is  skin 
friction.  The  skin  friction  on  an  object  is  a  function  of  the  fluid  density  p,  velocity  U,  drag 
coefficient  Cd,  and  the  surface  area  A. 

Fd  =  ~pU2CDA  (3.14) 
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The  drag  coefficient  is  a  function  of  the  Reynolds  number.  The  Reynolds  number 
is  a  function  of  the  fluid  density,  velocity,  the  length  over  which  the  flow  varies  (the 
width  of  the  vane  in  this  case)  W,  and  the  dynamic  viscosity  of  the  fluid  (p). 

pUW 


Re  = 


M 


(3.15) 


Stokes’  drag  occurs  at  very  low  Reynolds  numbers  (Re  <  1).  The  Reynolds 
numbers  for  every  trial  in  the  horizontal  vane  radiometer  experiment  was  less  than  0.02. 
Therefore,  all  the  trials  in  this  experiment  were  securely  in  the  Stokes’  drag  regime.  The 
drag  coefficient  for  Stokes’  drag  is  proportional  to  the  inverse  of  the  Reynolds  number 
where  8 is  a  numerical  factor. 


CD  = 


SjLl 

pUW 


(3.16) 


For  thin  plates  oriented  parallel  to  the  flow,  a  theoretical  value  of  4.12  for  8  has 
been  calculated  [18].  On  the  other  hand,  an  experimental  value  of  5.91  has  been  observed 
[19].  Thus,  for  uniform  flow  past  a  flat  plate  the  Stokes  drag  force  on  both  sides  of  the 
plate  is 

Fd  =  2^pU2CDLW  =  5.91  UpL, 

where  L  is  the  dimension  perpendicular  to  the  flow  (the  length  of  the  vane  in  this  case). 


(3.17) 


The  Stokes’  drag  equation  does  not  take  edge  effects  into  account.  Each  edge 
parallel  to  the  flow  that  is  immersed  in  the  fluid  adds  the  empirically  determined 
correction  to  the  drag  force  equation  [20] 

Fedge=\.6UpW.  (3.18) 


At  low  pressure  there  is  an  apparent  decrease  in  the  dynamic  viscosity  due  to  the 
slip  length  [17] 


M  =■ 


pW 


W  +  2A 


2-a  ' 


a 


Chapman’s  expression  for  the  dynamic  viscosity  [16], 

p  =  —pCA, 

32 


(3.19) 


(3.20) 
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which  accounts  for  variable  collision  rate,  improves  on  Maxwell’s  classical  estimate, 
H  =  pCXf 3.  Here  C  is  the  mean  speed  of  a  gas  molecule.  Note  that  p  is  pressure 

independent.  The  substitutions  p=nm  and  C  =  yJSkHT /71m  ,  where  m  is  the  molecular 
mass,  express  //in  terms  of  kinetic  theory  quantities 


<rcsw  + 


yJmkBTW 
Z . V2 2-a 


(3.21) 


It  is  not  obvious  that  the  slip  length  correction  should  apply  to  the  drag  on  the  edges  but  a 
correction  is  necessary  for  the  edge  correction  to  the  drag  force  to  be  equal  to  zero  at  zero 
pressure.  We  will  assume  the  correction  takes  a  similar  form  for  both  the  skin  friction  and 
the  edge  drag. 


For  a  uniformly  rotating  vane,  including  edge  effects,  the  drag  force  is 

f 

Fd  =  co  Ls+-  //(5.91L  +  3.2W),  (3.22) 

V  2  J 

where  Ls  is  the  length  of  the  stem.  Figure  7  shows  the  drag  force  predicted  by  this 
equation  for  the  vanes  in  the  horizontal  vane  radiometer  experiment  at  10  rad/s. 


Figure  7.  Drag  Force  at  10  rad/s  (Theory) 


The  narrow  vanes  are  8  mm  by  16  mm  and  the  wide  vanes  are  16  mm  by  16  mm.  Both  vanes 
have  a  4  mm  stem  connecting  them  to  the  spindle  axis.  The  accommodation  coefficient  is  0.6. 
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c. 


ANGULAR  SPEED 


Assuming  negligible  mechanical  friction  in  the  spindle,  the  equilibrium  angular 
speed  co  of  the  horizontal  vane  radiometer  occurs  when  the  drag  torque  equals  the 
thermal  creep  torque.  Because  the  drag  force  is  linearly  related  to  the  angular  speed,  we 
can  define  the  drag  force  per  unit  angular  velocity  fD  as 

fr)=—-  (3-23) 

co 


Thus,  in  terms  of  thermal  creep, 


co  = 


T  F 

' tc1  TC 

rdfD 


(3.24) 


where  rtc  and  are  the  effective  moment  arms  of  the  thermal  creep  and  drag  forces 
respectively.  Please  note  that  in  our  experiment’s  geometry,  the  drag  force  moment  arm 
is  121%  and  125%  of  the  thermal  creep  force  moment  arm  for  the  narrow  and  wide  vane 
respectively  because  the  thermal  creep  force  acts  uniformly  along  the  length  of  the  vane 
while  the  drag  force  increases  further  from  the  rotation  axis.  Since  both  forces  will  scale 
linearly  with  the  number  of  vanes  it  is  possible  to  examine  the  torques  on  a  single  vane  to 
predict  the  angular  speed  of  a  radiometer  with  four  vanes.  Figure  8  shows  the  predicted 
angular  velocity  that  results  by  using  the  expressions  (3.13)  and  (3.22)  in  the  above 
equation. 
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Figure  8.  Angular  Speed  (theory) 
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The  narrow  vanes  are  8  mm  by  16  mm  and  the  wide  vanes  are  16  mm  by  16  mm.  Both 
vanes  have  a  4  mm  stem  connecting  them  to  the  spindle  axis.  The  temperature  difference 
on  both  vanes  is  9  K.  The  characteristic  length  of  the  temperature  gradient  on  the  surface 
is  3.5  mm.  The  accommodation  coefficient  is  0.6.  Positive  angular  speed  represents 
rotation  with  the  cold  side  leading.  The  angular  speeds  are  similar  at  10  mTorr. 
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IV.  HORIZONTAL  VANE  RADIOMETER  SIMULATIONS 


Direct  Simulation  Monte  Carlo  (DSMC)  is  a  method  of  determining  continuum 
properties  of  a  fluid  by  simulating  individual  particles  of  the  fluid  [21].  Simulation 
particles  are  weighted  to  represent  multiple  fluid  particles  to  reduce  computational  effort. 
The  Navier-Stokes  equations  of  continuum  gas  dynamics  are  a  computationally  easier 
method  of  determining  the  same  quantities  but  provide  inaccurate  results  when  the  mean 
free  path  of  the  gas  approaches  a  representative  length  scale  of  the  system.  For  this  reason 
the  DSMC  method  is  most  usefully  applied  to  rarefied  gases  where  the  mean  free  path  of 
the  gas  is  of  the  same  order  or  larger  than  a  characteristic  length  of  the  system. 

DSMC  programs  divide  the  simulation  space  into  grid  cells.  Particles  are  initially 
created  in  random  locations  with  a  Gaussian  velocity  distribution  centered  on  an  assigned 
fluid  velocity.  At  each  time  step,  particles  move  a  distance  equal  to  their  velocity  times 
the  time  step.  The  behavior  of  particles  that  reach  the  boundary  of  the  simulation  space  is 
determined  by  the  boundary  conditions  of  the  simulation.  Particles  that  impact  a  surface 
in  the  simulation  space  undergo  either  specular  or  diffuse  reflection  depending  on  the 
simulation  inputs.  Momentum  transfer  between  the  surface  and  colliding  particles  is 
recorded  for  subsequent  calculation  of  the  pressure  on  the  surface  with  the  mass  of  a 
particle  equaling  the  mass  of  a  molecule  times  the  number  of  molecules  represented  by  a 
particle  in  the  simulation. 

After  the  particles  have  completed  their  motion,  particles  are  randomly  selected 
for  collision.  The  probability  that  an  individual  particle  will  be  selected  for  collision  is 
determined  by  the  mean  free  time  between  collisions  of  molecules  in  the  gas  being 
simulated  divided  by  the  time  step.  Collision  partners  are  then  chosen  from  within  the 
same  grid  cell.  The  collision  partners  are  assigned  new  velocities  that  conserve  energy 
and  momentum. 

Statistical  fluctuations  in  DSMC  results  can  be  averaged  out  by  decreasing  the 
number  of  molecules  represented  per  simulation  particle,  increasing  the  amount  of  time 
simulated,  or  averaging  the  results  of  multiple  trials.  Systemic  errors  can  be  introduced 
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by  large  grid  sizes  or  time  steps.  The  truncation  error  introduced  by  the  grid  size  is 
proportional  to  the  ratio  squared  of  the  length  of  a  side  of  the  grid  cell  to  the  mean  free 
path  of  the  gas  [22].  The  truncation  error  introduced  by  the  time  step  is  proportional  to 
the  ratio  squared  of  the  time  step  to  the  mean  free  time  between  collisions  of  molecules  in 
the  gas  [23]. 

SPARTA  (Stochastic  PArallel  Rarefied-gas  Time-accurate  Analyzer)  is  an  open- 
source  DSMC  program  developed  at  Sandia  National  Laboratories  [24].  We  used  it  in 
this  research  to  determine  separately  the  shear  pressures  caused  by  drag  and  thermal 
creep  on  vanes  similar  to  those  in  the  horizontal  vane  radiometer  experiment  at  the  range 
of  pressures  tested  in  the  experiment.  We  used  the  February  21,  2015  version  of 
SPARTA  with  one  modification.  The  modification  was  the  removal  of  the  command  to 
terminate  the  trial  if  a  particle  collides  with  an  interior  surface  of  the  vane.  This  feature 
was  in  the  code  to  ensure  objects  in  the  simulation  space  are  airtight,  but  occasionally  the 
algorithm  used  by  SPARTA  to  distinguish  between  interior  and  exterior  collisions 
mistakenly  classifies  a  collision  in  our  geometry.  We  checked  the  simulated  vanes  for 
airtightness  before  running  the  trials  on  the  modified  program.  We  ran  the  trials  on  the 
Hamming  supercomputer  at  the  Naval  Postgraduate  School  and  the  Spirit  supercomputer 
at  the  Air  Force  Research  Laboratory.  Each  trial  took  five  thousand  to  50  thousand  CPU 
hours  to  complete. 

A.  THERMAL  CREEP  SIMULATIONS 

We  ran  a  total  of  30  thermal  creep  simulations  encompassing  both  the  narrow 
vane  and  the  wide  vane  at  the  15  pressures  tested  in  the  experiment.  In  each  trial,  we 
simulated  a  single  stationary  vane  and  the  shear  pressure  on  the  vane  was  the  output.  We 
divided  each  trial  into  10  time  segments  which  reported  10  separate  average  shear 
pressures.  This  confirmed  the  behavior  was  not  transitory  and  allowed  us  to  calculate  a 
standard  deviation  between  the  10  segments  to  quantify  the  uncertainty  in  the  output. 

In  SPARTA,  geometry  objects  have  a  single  assigned  temperature  and  cannot 
touch  other  objects.  To  simulate  a  vane  with  a  temperature  gradient,  we  created  nine 
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objects  in  the  simulation  space.  Figure  9  is  an  overhead  view  of  the  wide  and  narrow 
vane  as  modeled  in  SPARTA. 


Figure  9.  Vanes  as  Modeled  in  SPARTA 


The  wide  and  narrow  vanes  are  each  divided  into  nine  segments  separated  by  two 

Angstroms.  Each  segment  has  a  single  assigned  temperature. 

The  cold  side  of  the  vanes  was  assigned  a  temperature  of  293  K.  The  hot  side  of  the 
vanes  was  assigned  a  temperature  of  302  K.  Each  intermediate  segment  was  1.125  K 
warmer  than  its  colder  neighbor.  All  nine  rectangular  boxes  which  comprised  a  vane  were 
16  mm  long  and  0.1  mm  thick.  The  seven  boxes  in  the  gradient  region  were  0.5  mm  wide 
which  from  equation  (3.10)  is  thin  enough  to  create  a  smooth  gradient  in  the  gas  at  all 
pressures  simulated.  The  two  end  boxes  were  6.25  mm  wide  for  the  wide  vane  and  2.25 
mm  wide  for  the  narrow  vane. 

The  simulation  space  was  18  mm  by  26  mm  by  10  mm  for  the  narrow  vane 
simulations  and  26  mm  by  26  mm  by  10  mm  for  the  wide  vane  simulations.  This  ensured 
at  least  one  mean  free  path  between  the  vane  and  the  boundary  at  all  pressures  simulated. 
The  boundary  condition  was  inflow  and  outflow  so  particles  whose  trajectory  took  them 
outside  the  simulation  space  would  no  longer  be  simulated  and  new  particles  would  be 
created  and  entered  into  the  simulation  space  throughout  the  trial.  The  size  of  the  grid 
cells  was  varied  with  pressure  so  that  each  edge  of  a  grid  cell  equaled  half  of  a  mean  free 
path  in  length.  The  total  number  of  grid  cells  ranged  from  280  for  the  narrow  vane  at  10 
mTorr  to  over  6  million  for  the  wide  vane  at  250  mTorr. 

Particles  were  created  with  an  initial  average  temperature  of  293  K.  Eight- 

hundred  and  fifty  million  particles  were  initially  created  in  the  narrow  vane  simulations 
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and  1.2  billion  particles  were  initially  created  in  the  wide  vane  simulations.  This  ensured 
there  was  an  average  of  200  particles  per  grid  cell  at  the  highest  pressure  trials.  The  lower 
pressure  trails  had  more  than  200  particles  per  grid  cell  since  the  grid  cells  were  larger 
and  fewer  in  the  lower  pressure  trials.  The  number  of  molecules  represented  by  a  single 
simulation  particle  varied  with  pressure.  This  value  varied  from  1.8  million  molecules  per 
simulation  particle  at  10  mTorr  to  45  million  molecules  per  simulation  particle  at  250 
mTorr.  79%  of  the  particles  were  diatomic  Nitrogen  and  21%  of  the  particles  were 
diatomic  Oxygen.  We  used  the  variable  soft  sphere  collision  model. 

The  time  step  was  100  ns.  This  value  is  one  fourth  the  mean  free  time  between 
collisions  at  250  mTorr  and  a  smaller  fraction  of  the  mean  free  time  between  collisions  at 
lower  pressures.  Each  simulation  ran  100  thousand  time  steps.  We  calculated  the  force  by 
multiplying  the  shear  pressure  output  by  the  area  of  the  surfaces.  Figure  10  shows  the 
results  of  this  calculation. 


Figure  10.  Thermal  Creep  Force  (Simulation) 
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The  narrow  vane  is  8  mm  by  16  mm  and  the  wide  vane  is  16  mm  by  16  mm.  The 
temperature  difference  on  both  vanes  is  9  K.  The  characteristic  length  of  the  temperature 
gradient  on  the  surface  is  3.5  mm.  The  accommodation  coefficient  is  0.6.  The  range  bars 
are  the  length  of  the  standard  deviation  of  the  10  segments  of  the  trial  at  the  reported 
pressure.  Range  bars  that  are  not  visible  are  hidden  behind  the  identifying  symbol.  A 
positive  force  is  directed  from  the  hot  side  to  the  cold  side.  We  believe  the  results  at  10 
mTorr  to  be  spurious. 
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It  is  interesting  to  note  that  the  values  are  negative  at  10  mTorr  (negative  meaning 
pointing  toward  the  hot  side  so  that  the  radiometer  would  rotate  with  the  black  side 
leading).  For  the  wide  vane,  the  magnitude  of  the  force  is  less  than  the  standard  deviation. 
For  the  narrow  vane  though,  the  force  is  -2.57  nN  with  a  1.63  nN  standard  deviation. 
None  of  the  10  time  segments  of  the  trial  reported  a  positive  force.  This  result  is  contrary 
to  the  experiment  and  theory,  but  similar  results  have  been  found  by  other  researchers 
using  the  DSMC  method  to  find  thermal  creep  forces  [1],  [2].  To  test  the  validity  of  this 
result  we  reran  the  simulation  at  10  mTorr  with  the  narrow  vane  in  the  larger  simulation 
space  used  previously  for  the  wide  vane.  The  reported  force  was  2.31  nN  in  the  positive 
direction  (pointing  toward  the  cold  side).  This  demonstrates  the  result  is  sensitive  to  the 
size  of  the  simulation  space  at  10  mTorr  where  there  was  a  single  mean  free  path  between 
the  vane  and  the  boundary  of  the  simulation  space  and  the  negative  results  in  our  trials 
are  not  to  be  trusted.  We  reran  the  simulation  in  the  same  manner  at  50  mTorr  where 
there  were  five  mean  free  paths  between  the  vane  and  the  boundary  of  the  simulation 
space.  The  reported  force  was  36.61  nN  in  the  positive  direction  which  compares 
favorably  with  the  previously  determined  force  of  35.55  nN  in  the  positive  direction  with 
a  9.64  nN  standard  deviation.  This  suggests  the  results  are  stable  with  respect  to  the  size 
of  the  simulation  space  when  there  are  multiple  mean  free  paths  between  the  vane  and  the 
boundary  of  the  simulation  space.  The  trials  reported  here  at  20  mTorr  and  above  include 
more  than  one  mean  free  path  between  the  vane  and  the  boundary  of  the  simulation 
space. 

We  also  ran  trials  at  selected  pressures  with  a  higher  accommodation  coefficient. 
A  higher  accommodation  coefficient  resulted  in  larger  force  values  as  expected.  Figure 
1 1  shows  these  results. 
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Figure  11.  Thermal  Creep  Force  with  Varied  Accommodation  (Simulation) 
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The  narrow  vanes  are  8  mm  by  16  mm  and  the  wide  vanes  are  16  mm  by  16  mm.  The 
temperature  difference  on  both  vanes  is  9  K.  The  characteristic  length  of  the  temperature 
gradient  on  the  surface  is  3.5  mm. 

B.  DRAG  FORCE  SIMULATIONS 

We  ran  drag  force  simulations  identical  to  the  thermal  creep  simulations  with  two 
exceptions.  First,  we  modeled  the  vanes  as  a  single  rectangular  object  with  a  temperature 
of  293  K.  Second,  collisions  between  the  vane  and  particles  were  calculated  as  if  the  vane 
were  rotating  at  10  radians  per  second  around  an  axis  4  mm  from  the  edge  of  the  vane 
though  the  vane  was  not  actually  rotating  in  the  simulation  space.  Figures  12  and  13  show 
the  results  of  these  trials. 
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Figure  12.  Drag  Force  at  10  rad/s  (Simulation) 


The  narrow  vane  is  8  mm  by  16  mm  and  the  wide  vane  is  16  mm  by  16  mm.  Particle- 
surface  collisions  are  calculated  as  if  the  vane  is  rotating  about  a  point  4  mm  beyond  its 
edge.  The  accommodation  coefficient  is  0.6. 


Figure  13.  Drag  Force  with  Varied  Accommodation  (Simulation) 
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The  narrow  vane  is  8  mm  by  16  mm  and  the  wide  vane  is  16  mm  by  16  mm.  Particle- 
surface  collisions  are  calculated  as  if  the  vane  is  rotating  at  10  rad/s  about  a  point  4  mm 
beyond  its  edge. 
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Additionally,  we  ran  a  number  of  drag  force  simulations  with  the  segmented 
vanes  used  in  the  thermal  creep  simulations  to  observe  the  effect  of  the  segmentation.  As 
seen  in  Figure  14,  the  segmentation  does  not  appreciably  change  the  drag  force  results 
leading  us  to  assume  the  segmentation  did  not  appreciably  change  the  thermal  creep  force 
results  either. 


Figure  14.  Drag  Force  with  Segmented  and  Unsegmented  Vanes  (Simulation) 
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The  effect  of  segmentation  is  negligible  at  10  and  50  mTorr.  In  our  simulation  setup, 
statistical  fluctuations  are  greater  at  higher  pressures  explaining  the  slight  differences  at 
100  mTorr. 
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c. 


ANGULAR  SPEED 


We  calculated  the  angular  speed  predicted  by  the  simulations  in  the  same  manner 
as  we  calculated  the  angular  speed  predicted  by  theory, 


co  =  co0 


(4.1) 


where  coo  is  the  angular  speed  in  the  simulation.  This  method  is  possible  because  the  drag 
force  is  linearly  proportional  to  velocity  at  these  pressures  as  seen  in  equation  (3.22)  and 
both  forces  scale  linearly  with  the  number  of  vanes  on  the  radiometer.  Figures  15  and  16 


report  the  angular  speed  predicted  by  the  simulations.  We  omitted  the  results  at  10  mTorr 
because  of  the  questionable  nature  of  the  results  of  our  thermal  creep  simulations  at  that 


pressure. 


Figure  15.  Angular  Speed  (Simulation) 


The  narrow  vane  is  8mm  by  16mm  and  the  wide  vane  is  16  mm  by  16  mm.  Both  vanes 
have  a  4  mm  stem  connecting  them  to  the  spindle  axis.  The  temperature  difference  on 
both  vanes  is  9  K.  The  characteristic  length  of  the  temperature  gradient  on  the  surface  is 
3.5  mm.  The  accommodation  coefficient  is  0.6.  Positive  angular  speed  represents  rotation 
with  the  cold  side  leading. 
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Figure  16.  Angular  Speed  with  Varied  Accommodation  (Simulation) 
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The  narrow  vane  is  8  mm  by  16  mm  and  the  wide  vane  is  16  mm  by  16  mm.  Both  vanes 
have  a  4  mm  stem  connecting  them  to  the  spindle  axis.  The  temperature  difference  on 
both  vanes  is  9  K.  The  characteristic  length  of  the  temperature  gradient  on  the  surface  is 
3.5  mm.  Positive  angular  speed  represents  rotation  with  the  cold  side  leading. 
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V.  COMPARISON  OF  EXPERIMENT,  THEORY,  AND 

SIMULATIONS 


In  both  the  simulations  and  theory  of  the  thermal  creep  force,  the  force  rises  as  a 
function  of  pressure  at  very  low  pressures  then  becomes  almost  constant  as  a  function  of 
pressure  at  relatively  higher  pressures.  Quantitatively,  the  differences  at  most  pressures 
are  less  than  an  order  of  magnitude.  A  qualitative  difference  does  stand  out,  though.  The 
simulations  show  the  force  on  the  narrow  vane  and  the  wide  vane  asymptotically 
approaching  different  values  as  the  pressure  increases  while  the  theory  has  the  force  on 
the  vanes  asymptotically  approaching  similar  values.  Figure  17  compares  the  simulation 
results  and  theory  predictions  at  the  range  of  pressures  examined.  We  omitted  the 
simulation  data  points  at  10  mTorr  because  of  the  questionable  validity  of  our  simulation 
parameters  at  that  pressure. 


Figure  17.  Thermal  Creep  Force  (Simulation  and  Theory) 
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The  narrow  vanes  are  8  mm  by  16  mm  and  the  wide  vanes  are  16  mm  by  16  mm.  The 
temperature  difference  on  the  vanes  is  9  K.  The  characteristic  length  of  the  temperature 
gradient  on  the  surface  is  3.5  mm.  The  accommodation  coefficient  is  0.6.  Positive  force  is 
directed  from  the  hot  side  to  the  cold  side. 
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In  both  the  simulations  and  the  theory  of  the  drag  force,  the  force  rises  with 
pressure  but  the  rate  of  increase  decreases  with  pressure.  The  change  in  the  rate  of 
increase  is  more  pronounced  in  the  theory  than  the  simulations.  Quantitatively  the 
agreement  is  very  good  around  100  mTorr  and  the  differences  are  less  than  a  factor  of 
two  at  all  pressures  examined.  Figure  18  compares  the  simulation  results  and  theory 
predictions  at  the  range  of  pressures  examined. 


Figure  18.  Drag  Force  (Simulation  and  Theory) 
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The  narrow  vanes  are  8  mm  by  16  mm  and  the  wide  vanes  are  16  mm  by  16  mm.  The 
vanes  have  a  4  mm  stem  connecting  them  to  the  spindle  axis.  They  are  rotating  at  10 
rad/s.  The  accommodation  coefficient  is  0.6. 


The  experimental  results,  simulation  predictions,  and  theory  predictions  of  the 
angular  speed  of  both  radiometers  have  a  peak  as  a  function  of  pressure  between  20 
mTorr  and  90  mTorr.  Quantitatively  the  differences  at  most  pressures  are  less  than  an 
order  of  magnitude.  The  theory  is  better  than  the  simulations  at  predicting  the 
experimental  results  for  the  narrow  vane  and  the  opposite  is  true  for  the  wide  vane.  The 
decrease  in  the  angular  speed  at  relatively  higher  pressures  is  more  pronounced  in  the 
experiment  than  in  the  simulation  or  theory  predictions.  The  difference  between  the 
asymptotic  force  values  on  the  narrow  vane  and  the  wide  vane  in  the  thermal  creep 
simulations  causes  the  simulation  prediction  of  the  angular  speed  to  be  greater  for  the 
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wide  vane  than  the  narrow  vane.  The  theory  prediction  has  the  narrow  vane  faster  than 
the  wide  vane  at  most  pressures.  The  experimental  results  qualitatively  match  the  theory 
predictions  in  this  regard.  Figures  19  and  20  compare  the  experimental  results  and 
simulation  and  theory  predictions  at  the  range  of  pressures  examined. 


Figure  19.  Narrow  Vane  Angular  Speed  (Experiment,  Simulation,  and  Theory) 
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The  vanes  are  8  mm  by  16  mm.  The  radiometer  has  4  vanes  attached  to  the  spindle  by  4 
mm  stems.  The  temperature  difference  on  the  vanes  is  9  K.  The  characteristic  length  of 
the  temperature  gradient  on  the  surfaces  is  3.5  mm.  Positive  angular  speed  represents 
rotation  with  the  cold  side  leading. 
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Figure  20.  Wide  Vane  Angular  Speed  (Experiment,  Simulation,  and  Theory) 
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The  vanes  are  16  mm  by  16  mm.  The  radiometer  has  4  vanes  attached  to  the  spindle  by  4 
mm  stems.  The  temperature  difference  on  the  vanes  is  9  K.  The  characteristic  length  of 
the  temperature  gradient  on  the  vanes  is  3.5  mm.  Positive  angular  speed  represents 
rotation  with  the  cold  side  leading. 


The  general  qualitative  and  better  than  order  of  magnitude  quantitative  agreement 
among  the  results  of  the  three  methods  of  investigation  provides  confidence  that  each 
method  of  investigation  is  sound.  Typically,  DSMC  results  show  better  quantitative 
agreement  with  experimental  results  when  the  simulation  design  and  parameters  more 
closely  resemble  the  actual  conditions  of  the  experiment,  but  this  is  often  at  the  expense 
of  greater  computational  resources.  Factors  which  tend  to  cause  our  simulations  to  show  a 
certain  degree  of  quantitative  disagreement  with  the  experimental  results  include  the 
method  used  to  model  the  gradient,  running  separate  thermal  creep  and  drag  force 
simulations,  running  simulations  on  static  vanes  thereby  negating  drafting  and  other 
aerodynamic  effects,  and  simulating  only  the  immediate  volume  of  air  surrounding  the 
vane  rather  than  the  entire  chamber.  Eliminating  these  factors  would  have  made  the 
simulations  time  and  cost  excessive  if  not  unfeasible.  In  our  theoretical  results,  we  chose 
to  use  reasonable  assumptions  for  the  numerical  factors  a,  /?,  and  x,  but  these  values  could 
be  adjusted  to  better  fit  the  experimental  and  simulation  data.  For  example,  Figure  21 
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shows  the  agreement  between  the  narrow  vane  results  if  a  value  of  0.33  is  assumed  for  [> 
and  .r. 


Figure  21.  Narrow  Vane  Angular  Speed  (ft  =  x  =  0.33) 
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The  vanes  are  8  mm  by  16  mm.  The  radiometer  has  4  vanes  attached  to  the  spindle  by  4 
mm  stems.  The  temperature  difference  on  the  vanes  is  9  K.  The  characteristic  length  of 
the  temperature  gradient  on  the  surfaces  is  3.5  mm.  Positive  angular  speed  represents 
rotation  with  the  cold  side  leading.  The  values  of  certain  numerical  factors  in 
the  theoretical  results  were  adjusted  to  more  closely  match  the  experimental  data. 
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VI.  MICRON-SCALE  DEVICES  AT  ATMOSPHERIC  PRESSURE 


Microelectromechanical  System  (MEMS)  devices  can  experience  rarified  gas 
effects  at  atmospheric  pressure  as  their  dimensions  approach  the  mean  free  path  of  air  at 
atmospheric  pressure  (68  nm).  In  the  experiment,  we  observed  thermal  creep  forces  even 
when  the  width  of  the  vane  was  80  times  larger  than  the  mean  free  path  of  the 
surrounding  air.  This  presents  possibilities  for  applications  of  the  thermal  creep  force  on 
micron-scale  devices  at  atmospheric  pressure. 

A.  SCALING 

As  seen  in  equation  (3.13)  the  thermal  creep  force  scales  linearly  with  the 
temperature  difference  on  the  surface.  It  is  therefore  necessary  to  understand  the  scaling 
relationship  between  the  temperature  differences  on  a  millimeter-scale  vane  and  a 
micron-scale  vane. 

Before  a  light  is  shone  on  a  radiometer  vane,  the  two  halves  are  at  thermal 
equilibrium  with  the  environment  and  each  other.  The  apparent  temperature  of  the 
environment  increases  when  the  vane  is  exposed  to  a  new  light  source.  Both  halves  of  the 
vane  will  increase  in  temperature  until  they  are  again  at  equilibrium  with  the  environment 
and  each  other,  but  the  black  side  of  the  vane  will  initially  increase  in  temperature  more 
quickly.  As  a  function  of  time,  then,  the  temperature  difference  between  the  two  halves  of 
the  vane  will  start  at  zero  before  being  exposed  to  a  new  light  source,  will  increase  to  a 
maximum  as  the  black  side  heats  up  more  quickly,  will  decrease  as  the  white  side  catches 
up  to  the  black  side,  and  will  then  again  be  zero  as  both  halves  are  again  in  equilibrium 
with  the  environment. 

The  temperature  difference  on  a  vane  AT  is  a  function  of  the  thermal  energy  of  the 
two  halves  where  Th  is  the  temperature  of  the  hot  side,  7c  is  the  temperature  of  the  cold 
side,  C  is  the  heat  capacity  assumed  to  be  the  same  on  both  halves,  EH  is  the  thermal 
energy  of  the  hot  side,  and  Ec  is  the  thermal  energy  of  the  cold  side. 

AT=Th-Tc=±(Eh-Ec)  (6,1) 
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The  temperature  difference  is  a  maximum  when  the  time  rates  of  change  of  the 
thermal  energies  of  the  two  halves  are  equal. 


a(A7~)  0  ,beh  sec 

dt  dt  dt 


(6.2) 


The  temperature  difference  on  the  vanes  in  the  experiment  was  not  observed  to  vary  with 
pressure  suggesting  convective  heat  transfer  between  the  two  sides  is  negligible.  The  time 
rates  of  change  of  the  thermal  energies  are  thus  functions  of  their  respective  net  black 
body  radiation  and  the  thermal  conduction  rate  between  them  where  cr5B  is  the  Stephan- 
Boltzmann  constant,  s  is  the  emissivity,  Aexp  is  the  surface  area  of  one  half  exposed  to  the 
new  light  source,  Tenv  is  the  apparent  temperature  of  the  environment,  k,  is  the  thermal 
conductivity  between  the  two  halves,  Acon  is  the  area  in  contact  between  the  two  halves  of 
the  vane,  and  Wo  is  the  width  of  the  temperature  gradient  in  the  vane. 


dE^ 

dt 

dEc 

dt 


AT 

Wo 

AT 


=  Vc\,V(rL-Tc)  +  K\on  — 


(6.3) 

(6.4) 


The  temperature  of  the  hot  side  of  the  vane  can  be  rewritten  as  a  function  of  the 
temperature  of  the  cold  side  and  the  temperature  difference. 


th=tc+at=tc 


1+ 


AT 


lc  J 


(6.5) 


For  the  vanes  in  the  experiment  and  most  applications,  the  temperature  difference  is 
much  smaller  than  the  temperature  of  the  cold  side  of  the  vane  allowing  the  following 
approximation  to  be  made. 


Ti  =  T4 

H  C 


1  + 


AT 


[c  J 


:T4+at(;at 


(6.6) 


Substituting  approximation  (6.6)  into  equation  (6.3)  then  setting  equal  to  equation 
(6.4)  and  solving  for  the  temperature  difference  provides  the  following  result. 


AT 
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AfHzi MJtLzE) 

4(7  s  A  T3  +  A 

SBb H  ^exp1  C  ^  it, 

wG 


(6.7) 
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The  scaling  behavior  of  this  relation  depends  on  whether  the  blackbody  term, 
4(7 Sf K  jl >  or  the  thermal  conductivity  term,  2 ktAcon/WG  ,  dominates  the  denominator. 


For  the  vanes  in  the  experiment,  the  blackbody  term  is  more  than  10  times  greater  than 
the  thermal  conductivity  term  due  to  their  relatively  large  areas  exposed  and 
characteristic  lengths  of  their  temperature  gradients  and  relatively  small  thermal 
conductivity  and  areas  of  contact  between  the  two  halves. 


mm-scale,  low  kt  vanes:  A T  « 


(f^2f c)KcIe1 

*eBTl 


(6.8) 


This  expression  also  confirms  that  the  magnitude  of  the  temperature  difference  was  not  a 
function  of  the  width  of  the  vanes  in  the  experiment. 


The  blackbody  term  is  a  function  of  two  factors  of  physical  size  ( Aexp  or  length 
times  width)  while  the  thermal  conductivity  term  is  a  function  of  one  factor  of  physical 
size  and  an  aspect  ratio  (Aco,/WG  or  length  times  thickness  divided  by  width).  For  micron- 
scale  vanes  then  the  thermal  conductivity  term  is  likely  to  dominate  the  denominator.  The 
smaller  dimensions  will  also  result  in  the  characteristic  length  of  the  temperature  gradient 
being  equal  to  the  width  of  the  vane.  Assuming  the  top  and  bottom  of  the  vane  are  both 
exposed  to  the  light  source,  the  following  is  an  approximate  expression  for  the  maximum 
temperature  difference  on  a  MEMS  vane  where  W  is  the  width  of  the  vane  and  H  is  the 
height  or  thickness  of  the  vane. 


MEMS  vanes:  A Tmax 


2 k,H 


(6.9) 


or 


AT . oc 


W~ 
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or 


AT 


■  oc  ■ 


W 


(6.10) 


(6.11) 


W  H 

The  temperature  gradient  is  thus  a  function  of  the  aspect  ratio  of  the  width  (which 
affects  the  net  blackbody  radiation)  to  the  thickness  (which  affects  the  thermal 
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conductivity).  It  is  not  a  function  of  the  length  (which  affects  both).  The  temperature 
difference  is  simply  the  roughly  constant  temperature  gradient  times  the  width. 

Ignoring  slip  length  and  edge  corrections  for  the  purpose  of  developing  a  scaling 
relation,  the  thermal  creep  force  from  equation  (3.13)  is  linearly  proportional  to  the 
length  and  temperature  difference  of  the  vane. 

Ftc  cc  ATL  (6.12) 

Combining  this  result  with  equation  (6.10)  provides  the  following  scaling 
relationship  for  the  thermal  creep  force  on  a  MEMS  vane. 

W2L 

Ftc*—  (6-13) 

LI 

These  scaling  relationships  are  useful  in  two  ways.  Firstly,  they  allow  results  from 
centimeter-scale  experiments  to  be  extended  into  the  MEMS  realm.  Secondly,  they 
indicate  the  effect  of  varying  the  aspect  ratio  of  a  device.  For  a  systematic  use  of  the 
scaling  relations,  consider  shrinking  the  horizontal  vanes,  while  maintaining  geometrical 
similarity  (constant  aspect  ratios  between  the  physical  dimensions  of  the  vanes)  and 
increasing  the  pressure  to  its  atmospheric  value  while  maintaining  a  constant  ratio 
between  the  physical  dimensions  of  the  vanes  and  the  mean  free  path.  For  the  scaling 
relations  to  be  meaningful,  the  important  parameters  to  determine  experimentally  are  the 
temperature  difference  across  the  vane,  the  characteristic  length  of  the  temperature 
gradient  on  the  vane,  and  the  value  at  which  thermal  creep  was  observed.  From  the 
pressure  reading,  the  physical  dimensions  of  the  vanes  can  be  defined  as  functions  of  the 
mean  free  path.  Then,  equation  (3.13)  can  be  used  to  calculate  the  force  at  the  pressure 
where  the  width  of  the  vane  equals  the  characteristic  length  of  the  temperature  gradient 
on  the  surface.  Scaling  relationships  (6.10)  and  (6.13)  can  then  be  used  to  extend  the 
result  of  that  calculation  to  estimate  the  temperature  difference  and  force  on  a  MEMS 
device  at  atmospheric  pressure.  Finally,  the  aspect  ratios  in  the  device  can  be  modified  to 
achieve  a  desired  force. 

To  illustrate  this  process,  consider  the  narrow  vanes  in  the  experiment.  We 
observed  a  9  K  temperature  difference,  the  temperature  gradient  had  a  characteristic 
length  of  3.5  mm,  and  thermal  creep  was  observed  at  pressures  around  225  mTorr.  Thus, 
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the  width  of  the  vane  at  this  pressure  is  35  mean  free  paths  and  the  length  is  70  mean  free 
paths.  The  width  of  the  vane  equals  the  characteristic  length  of  the  temperature  gradient 
when  the  vane  is  3.5  mm  by  7  mm  and  the  mean  free  path  is  0.1  mm.  Using  these  values 
in  equation  (3.13)  provides  a  force  of  173  nN.  Scaling  relationships  (6.10)  and  (6.13)  then 
predict  a  6  mK  temperature  difference  and  a  0.08  pN  force  on  a  2.4  pm  by  4.8  pm  device 
at  atmospheric  pressure.  If  we  instead  wanted  a  force  of  1  pN,  relation  (6.13)  shows  we 
could  increase  the  length  to  60  pm  or  increase  the  width  to  8.5  pm. 

Ignoring  the  slip  length  correction  to  the  dynamic  viscosity  for  the  purpose  of 
developing  a  drag  force  scaling  relation  and  assuming  a  fixed  aspect  ratio,  the  drag  force 
on  a  rotating  vane  from  equation  (3.22)  is  proportional  to  a  measure  of  length  squared. 

FdocL2  (6.14) 

Since  both  the  thermal  creep  and  drag  forces  scale  as  a  measure  of  length  squared, 
devices  of  different  sizes  should  be  able  to  achieve  comparable  angular  velocities  at  their 
peak  pressures  though  the  peak  pressure  will  scale  with  the  size  of  the  device. 

Substituting  the  angular  velocity  U  =  coL  and  Maxwell’s  estimate  of  the  dynamic 
viscosity  //  =  /?C/l/3  into  the  Reynolds  number  Re  =  pUW / //  provides  the  following 
scaling  relation  for  the  Reynold’ s  number  for  horizontal  vane  radiometers  operating  at  the 
same  angular  speed  in  the  same  gas  at  the  same  temperature, 

Reoc— ,  (6.15) 

A 

confirming  a  MEMS  horizontal  vane  radiometer  will  remain  in  the  Stokes’  drag  regime  at 
much  higher  pressure  than  a  millimeter  scale  device. 

B.  DESIGN 

Thermal  creep  forces  will  occur  in  MEMS  devices  whenever  a  temperature 
gradient  in  the  surrounding  gas  is  parallel  to  a  surface.  A  geometrically  symmetrical 
component  that  is  hotter  than  the  ambient  temperature  will  create  temperature  gradients 
in  the  surrounding  gas  and  experience  thermal  creep  forces,  but  they  will  tend  to  cancel 
due  to  the  symmetry.  Thermal  creep  forces  that  are  likely  to  degrade  device  performance 
are  created  by  unevenly  heated  components  or  hot  components  with  an  asymmetrical 
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geometry  and  should  be  considered  in  the  design  process.  If  properly  designed,  MEMS 
devices  could  even  utilize  thermal  creep  forces  for  energy  harvesting  and  for  sensing. 
Passian  et  al.  designed  a  MEMS  cantilever  with  pores  that  experiences  thermal 
transpiration  forces  when  exposed  to  a  modulated  laser  [25].  When  the  modulation 
frequency  is  on  the  order  of  1  kHz  or  greater  the  cantilever  does  not  reach  a  uniform 
temperature  creating  a  temperature  gradient  in  the  structure.  Here  we  consider  two 
designs  based  on  our  horizontal  vane  radiometer  which  establishes  a  temperature  gradient 
by  having  two  halves  of  differing  emissivity. 

The  thermal  conductivity  of  silicon  is  2,600  times  greater  than  the  thermal 
conductivity  of  the  paper  that  was  used  in  the  horizontal  vane  radiometer  experiments.  As 
seen  in  equation  (6.9)  the  temperature  difference  created  on  a  MEMS  vane  is  inversely 
proportional  to  the  thermal  conductivity  making  low  thermal  conductivity  polymers  and 
ceramics  better  design  choices  than  silicon  or  metals.  Two  common  polymers  used  in 
MEMS  fabrication  are  the  polyimide  Kapton  and  the  epoxide  SU-8.  Analysis  of  the 
temperature  difference  on  a  radiometer  vane  made  of  two  materials  is  complicated  by  the 
differing  heat  capacity  and  thermal  conductivity  on  the  two  sides  of  the  vane,  but  the 
values  for  Kapton  and  SU-8  are  close  enough  that  the  scaling  analysis  can  be  used  for 
rough  estimation.  The  thermal  conductivities  of  Kapton  and  SU-8  are  2.4  and  4  times 
larger  than  the  thermal  conductivity  of  paper  respectively  [26],  [27],  but  gaps  in  the 
contact  area  between  the  two  materials  because  of  manufacturing  tolerances  will  decrease 
the  thermal  transfer  between  the  two  sides  of  the  vane.  Tabulated  emissivity  values  differ 
between  sources  and  are  difficult  to  find  for  certain  materials.  Emissivity  values  for  white 
paper  range  from  0.68  to  0.93.  Values  for  carbon,  which  provides  the  coloring  in  Laserjet 
ink,  range  from  0.80  to  0.95  [28].  The  emissivity  of  Kapton  is  0.57  and  the  emissivity  of 
SU-8  is  0.95  [29],  [30].  These  values  suggest  the  emissivity  difference  between  Kapton 
and  SU-8  would  be  even  greater  than  the  difference  between  white  paper  and  black 
LaserJet  ink.  These  material  properties  suggest  that  a  MEMS  horizontal  vane  radiometer 
made  of  polyimide  and  SU-8  could  rotate  with  an  angular  velocity  of  the  same  order  of 
magnitude  as  the  radiometers  in  our  experiment. 
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We  designed  a  MEMS  horizontal  vane  radiometer,  as  shown  in  Figure  22. 


Figure  22.  MEMS  Horizontal  Vane  Radiometer  Cross  Section  and  Vane 

Structure  Overhead  View 


Hub 

b— 


Vane  Spindle  yane 
Substrate 

The  figure  on  left  is  a  cross  section  of  a  MEMS  horizontal  vane  radiometer.  The  figure  on 
the  right  is  an  overhead  view  of  the  vane  structure. 

The  vane  structure  consists  of  a  piece  of  Kapton  and  a  separate  piece  of  SU-8.  Each  piece 
consists  of  a  circular  hub  and  four  vanes.  Each  piece  was  given  its  own  hub  to  make  the 
structure  more  robust.  The  spindle  has  a  base  and  a  top  both  of  which  are  wider  than  the 
openings  in  the  hubs  while  the  spindle  itself  is  narrower  than  the  hub  openings.  This 
design  prevents  the  vane  structure  from  riding  on  the  substrate  or  separating  from  the 
spindle.  We  used  MEMS  Pro  L-Edit  to  plan  the  steps  necessary  to  realize  the  design.  The 
fabrication  process  requires  11  steps.  Figures  23,  24,  and  25  show  a  selection  of  the 
fabrication  steps. 
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Figure  23.  MEMS  Horizontal  Vane  Radiometer  Fabrication  Step  5 


The  black  structure  is  the  SU-8  portion  of  the  radiometer.  The  blue  is  a  sacrificial  layer  of 
silicon  dioxide.  The  base  of  the  spindle  is  below  the  sacrificial  layer.  The  next  step  is  to 
deposit  a  layer  of  Kapton. 


Figure  24.  MEMS  Horizontal  Vane  Radiometer  Fabrication  Step  7 


The  white  is  the  Kapton  portion  and  the  black  is  the  SU-8  portion  of  the  radiometer.  The 
blue  is  a  sacrificial  later  of  silicon  dioxide.  The  base  of  the  spindle  is  below  the  sacrificial 
layer.  The  next  step  is  to  deposit  another  sacrificial  layer. 
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Figure  25.  MEMS  Horizontal  Vane  Radiometer 


The  white  is  the  Kapton  portion  and  the  black  is  the  SU-8  portion  of  the  radiometer.  The 
light  grey  circle  is  the  spindle  top.  The  vane  structure  is  resting  on  the  spindle  base  rather 
than  the  substrate.  A  narrow  spindle  runs  through  the  hub. 


C.  APPLICATIONS 

There  are  numerous  potential  applications  of  the  thermal  creep  force  on  MEMS 
devices.  One  of  the  most  interesting  is  energy  harvesting.  Consider  an  electric  generator 
with  a  horizontal  vane  radiometer  rotor  (see  Figure  26). 


Figure  26.  Simplified  Electric  Generator  with  Horizontal  Vane  Radiometer 

Rotor 


The  vanes  have  a  conductive  wire  along  their  axis.  Rotation  in  the  magnetic  field  creates  an 
electromotive  force  in  the  wire.  The  circuit  closes  when  the  vane  sweeps  the  sliding  contact. 
This  is  a  simplified  diagram  of  a  DC  generator  to  facilitate  calculation  of  the  voltage  produced. 
A  real  rotor  could  be  seated  inside  a  cylinder  with  the  vanes  in  contact  with  the  casing  so  the 
circuit  is  always  closed. 
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The  voltage  produced  is  proportional  to  the  angular  velocity,  magnetic  field  B ,  and  the 


length  of  the  vanes  squared  found  by  integrating  the  magnetic  force,  v  x  B ,  over  the 
length  of  a  vane. 


V  = 


coBL 2 
2 


(6.16) 


A  radiometer  with  100  pm  vanes  operating  a  10  rad/s  in  a  1  T  magnetic  field  would 
generate  50  nV.  The  small  size  though  would  allow  numerous  such  devices  to  fit  on  a 
single  chip.  1,000  connected  in  series  would  generate  50  pV. 


Alternatively,  consider  a  bicolored  cantilever  beam  with  a  thin  piezoelectric  film 
at  its  base  (see  Figure  27). 


Figure  27.  Cantilever  Beam  with  Piezoelectric  Base 


The  piezoelectric  film  is  in  grey  at  the  base  of  the  cantilever  beam.  The  deflection  of  the  beam  will  create  a 
tensile  stress  in  piezoelectric  film  along  its  length  which  will  generate  a  potential  difference  orthogonal  to 
the  stress. 
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The  voltage  produced  in  the  piezoelectric  film  is  as  follows  where  dj /  is  the  appropriate 
piezoelectric  coefficient,  F  is  the  total  force  assumed  to  be  applied  at  the  free  end  of  the 
cantilever  beam,  L  is  the  length  of  the  cantilever  beam,  W  is  the  width  of  the  cantilever 
bean,  t  is  the  thickness  of  the  piezoelectric  film,  e  is  the  permittivity  of  the  piezoelectric 
film,  and  /  is  the  area  moment  of  inertia  of  a  cross  section  of  the  cantilever  beam  [31]. 


_  d3lFLWt 
2sl 


(6.17) 


Substituting  WH/12  for  the  area  moment  of  inertia  where  H  is  the  height  or  thickness  of 
the  cantilever  beam  and  multiplying  by  3/8  to  account  for  the  fact  the  thermal  creep  force 
will  be  uniformly  distributed  over  the  cantilever  beam  and  not  applied  solely  at  the  free 
end  provides  the  following  result. 

9  d,,FLt 


V  = 


AsW2H 


(6.18) 


Combining  equations  (6.13)  and  (6.18)  provides  the  following  as  the  dependence  of  the 
voltage  on  the  dimensions  of  the  cantilever  beam. 


V  oc- 


H~ 


(6.19) 


The  voltage  is  thus  a  function  of  the  ratio  of  the  length  to  the  thickness  of  the  cantilever 
beam  rather  than  the  overall  size  and  completely  independent  of  the  width  of  the 
cantilever  beam.  The  devices  can  therefore  be  as  small  as  MEMS  manufacturing 
technology  will  allow  and  connected  in  series  to  create  a  higher  overall  voltage. 
According  to  equation  (6.18),  a  100  pm  by  10  pm  by  lpm  cantilever  beam  able  to 
produce  1  pN  of  thermal  creep  force  with  a  1  pm  thick  zinc  oxide  film  at  is  base  would 
produce  150  nV.  1,000  such  devices  connected  on  a  single  chip  would  produce  150  pV. 
The  disadvantage  of  this  design  compared  to  the  rotor  design  is  the  need  for  a  periodic 
interrupter  of  the  light  to  maintain  a  continuous  alternating  current. 
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VII.  CONCLUSION 


We  have  investigated  the  thermal  creep  shear  force  through  experiment,  kinetic 
theory  calculations,  and  the  DSMC  simulations.  All  three  investigations  show  the  force 
acting  from  the  hot  to  the  cold  side  of  the  surface.  The  theory  and  the  simulations  both 
show  the  force  rising  with  pressure  and  asymptotically  approaching  a  constant  value 
though  the  theory  shows  the  force  on  both  vanes  approaching  the  same  value  while  the 
simulations  show  the  force  on  the  two  vanes  approaching  different  values.  The  theory  and 
the  simulations  both  show  the  decrease  in  the  radiometer’s  angular  speed  as  a  function  of 
pressure  at  relatively  higher  pressures  to  be  the  result  of  increasing  drag  force  rather  than 
decreasing  thermal  creep  force. 

Left  unexplained  by  our  investigation  is  why  the  theory  is  better  at  predicting  the 
experimental  results  for  the  narrow  vane  and  the  simulation  results  are  better  at  predicting 
the  experimental  results  for  the  wide  vane.  This  suggests  that  at  least  two  of  the  three 
methods  of  investigation,  while  good,  can  be  further  improved  upon.  There  are  potential 
areas  of  future  work  which  would  provide  additional  insight  in  all  three  methods  of 
investigation.  In  our  experiment,  we  measured  angular  speed  which  is  the  result  of  all 
forces  on  the  moving  radiometer.  Additional  insight  on  the  thermal  creep  force  in 
particular  could  be  achieved  by  static  force  measurements  on  a  fixed  vane.  In  our  theory, 
we  made  a  quasi-static  assumption  and  made  assumptions  about  the  temperature 
distribution  of  the  gas.  Adding  dynamic  effects  to  the  theory  and  experimentally  or  in 
simulation  directly  observing  the  temperature  distribution  in  the  gas  and  incorporating  the 
results  could  improve  our  theory.  In  our  simulations,  the  thermal  creep  force  and  drag 
force  were  determined  in  separate  trials.  Trials  with  both  a  rotating  vane  and  a 
temperature  gradient  could  reveal  any  second  order  effects  from  the  interaction  of  drag 
and  thermal  creep. 

We  also  provided  scaling  relations  between  millimeter  scale  and  MEMS 
horizontal  vane  radiometers,  a  design  for  a  MEMS  horizontal  vane  radiometer,  and  two 
potential  applications  of  the  thermal  creep  force  for  energy  harvesting  with  MEMS 
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devices.  Future  work  to  optimize  these  applications  would  improve  the  modest  voltage 
output  we  predicted  and  possibly  create  useful  devices. 
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APPENDIX  A.  EINSTEIN  EFFECT 


In  Einstein’s  development  of  his  effect,  he  simplified  the  behavior  of  the  gas  by 
assuming  all  gas  molecules  moved  in  only  one  axis  direction  at  a  time.  In  a  cubic  volume, 
1/6  of  the  molecules  leave  the  volume  through  any  one  face  of  the  cube.  The  number  of 
molecules  impacting  a  surface  of  area  A  from  one  direction  per  unit  time  is  then 


—  nAv 

6 


(7.1) 


where  n  is  the  number  density  and  v  is  the  average  speed  of  a  gas  molecule.  The 
difference  in  momentum  between  molecules  of  mass  m  coming  from  the  hot  side  and  the 
cold  side  of  a  surface  perpendicular  to  a  temperature  gradient  is 

mv+  —  mv  (7.2) 

where  v+  and  v_  are  the  speeds  of  molecules  from  the  hot  and  cold  sides  respectively.  In  a 
gas  with  no  net  mass  flow,  the  number  of  molecules  impacting  the  surface  per  unit  time 
from  each  direction  is  equal.  The  force  on  the  surface  is  then  the  number  of  molecules 
impacting  from  one  direction  per  unit  time  times  the  difference  in  momentum  between 
molecules  coming  from  opposing  sides. 

F Einstein  =~  nAmv(v+  -  v_  )  (7.3) 

o 


where  the  temperature  gradient  is  in  the  positive  direction.  Rewriting  v  as  the  average  of 
v+  and  v.,  we  have 


F Einstein  =  ~\nAm 

o 


'V  +v  A 


V 


K-v_) 


Substituting  3^bT/2  =  otv2/2  we  have 


- nAmlvl  -  v 2 

12  1  +  " 


(7.4) 


Einstein  =  ~nkB{T+  ~T-)A- 


(7.5) 


Assuming  the  temperature  of  molecules  arriving  from  the  hot  and  the  cold  sides  were 
each  established  one  mean  free  path  away, 


T+  =T  +  A  — 
dx 


(7.6) 


and 
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(7.7) 


T  =T-y l  — 
dx 

where  dT/dx  is  the  temperature  gradient.  Substituting  equations  (7.6)  and  (7.7)  into 
equation  (7.5)  we  have 


1  l?dT 

k'B/ 


TX.  , .  =  -—nkBA——A. 


2  *  dx 

Then  substituting  the  ideal  gas  equation  where  P  is  pressure  we  have 

1  PA  dT 


Einstein 


2  T  dx 


-A. 


(7.8) 


(7.9) 


Finally,  assuming  that  the  temperature  gradient  is  proportional  to  the  temperature 
difference  between  the  hot  and  the  cold  sides  of  the  vane  AT  divided  by  the  mean  free 
path,  the  area  is  a  mean  free  path  wide  strip  the  length  of  the  perimeter  Lp,  and  dropping 
the  numerical  factor  as  this  is  a  qualitative,  order  of  magnitude  expression,  we  have  [11] 

FBi„„m=-P^Lr-  (7.10) 

In  terms  of  kinetic  theory  quantities, 

ATX 

FEinstein* - (7.11) 

&CS 

Scandurra  et  al.  find  a  similar  relationship  using  the  Chapman-Enskog  method  but 
include  a  factor  of  2 -a  to  include  recoil  effects  when  the  molecules  are  not  fully 
accommodated  by  the  surface.  [4] 


It  is  also  possible  to  relate  the  force  to  the  heat  flow  through  an  open  area  of  the 
same  size  and  orientation  to  the  temperature  gradient.  The  difference  in  energy  between 
molecules  coming  from  the  hot  side  and  the  cold  side  is 

— mvl-  —  mv2.  (7.12) 

2  2” 

The  heat  flow  is  then  the  number  of  molecules  passing  through  from  one  direction  per 
unit  time  times  the  energy  difference  between  molecules  coming  from  opposite  directions 
or 


6  A  =  ——nAv 

6 


—  mv, - mv 


-—nAmv{v2+  -  v2 


(7.13) 
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where  (|>  is  the  heat  flux.  Comparing  equation  (7.4)  to  (7.13)  makes  clear  that  [11] 


(7.14) 
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APPENDIX  B.  EXAMPLE  THERMAL  CREEP  INPUT  SCRIPT 


seed 

7092015 

dimension 

3 

global 

gridcut  0.001 

boundary 

0  0  0 

create_box 

-0.009  0.009  -0.013  0.013  -0.005  0.005 

create_grid 

35  50  19  block  *  *  * 

global 

temp  293 

global 

nrho  1.65e21  fnum  9.07e6 

global 

gravity  9.81  0  0-1 

species 

air.  species  N2  02 

mixture 

air  N2  frac  .79 

mixture 

air  02  frac  .21 

read_surf 

1  data.cold_vane8mm 

surf_collide 

1  diffuse  293  .6 

surf_modify 

collide  1  1 

read_surf 

2  data.mid_vanel 

surf_collide 

2  diffuse  294.125  .6 

surf_modify 

collide  2  2 

read_surf 

3  data.mid_vane2 

surf_collide 

3  diffuse  295.25  .6 

surf_modify 

collide  3  3 
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read_surf 

4  data.mid_vane3 

surf_collide 

4  diffuse  296.375  .6 

surf_modify 

collide  4  4 

read_surf 

5  data.mid_vane4 

surf_collide 

5  diffuse  297.5  .6 

surf_modify 

collide  5  5 

read_surf 

6  data.mid_vane5 

surf_collide 

6  diffuse  298.625  .6 

surf_modify 

collide  6  6 

read_surf 

7  data.mid_vane6 

surf_collide 

7  diffuse  299.75  .6 

surf_modify 

collide  7  7 

read_surf 

8  data.mid_vane7 

surf_collide 

8  diffuse  300.875  .6 

surf_modify 

collide  8  8 

read_surf 

9  data.hot_vane8mm 

surf_collide 

9  diffuse  302  .6 

surf_modify 

collide  9  9 

collide 

vss  air  air.vss 

collide_modify 

vremax  1000  yes 

create_particles 

air  n  0 

fix 

in  inflow  air  all 

timestep 

le-7 
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compute 

1  surf  air  shx  shy  shz 

fix 

1  ave/surf  1  100000  100000  c_l 

dump 

1  surf  100000  dump.vane50tcbb  id  vlx  vly  viz  v2x  v2y 

v2z  v3x  v3y  v3z  f_l 

compute 

2  surf  air  shx  shy  shz 

fix 

2  ave/surf  1  10000  10000  c_2 

dump 

2  surf  10000  dump.vane50tenthtcbb  id  vlx  vly  viz  v2x 

v2y  v2z  v3x  v3y  v3z  f_2 

stats 

1000 

stats_style 

step 

run 

100000 
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APPENDIX  C.  EXAMPLE  DRAG  INPUT  SCRIPT 


seed 

dimension 

global 

boundary 

create_box 

create_grid 

global 

global 

global 

species 

mixture 

mixture 

read_surf 

surf_collide 

surf_modify 

collide 

collide_modify 

create_particles 

fix 

time  step 
compute 
fix 


7092015 

3 

gridcut  0.001 
ooo 

-0.009  0.009  -0.013  0.013  -0.005  0.005 
35  50  19  block  *  *  * 
temp  293 

nrho  1.65e21  fnum  9.07e6 

gravity  9.81  0  0-1 

air.  species  N2  02 

air  N2  frac  .79 

air  02  frac  .21 

1  data. vane 

1  diffuse  293  .6  rotate  0-.0 12000  10 

collide  1  1 

vss  air  air.vss 

vremax  1000  yes 

air  n  0 

in  inflow  air  all 
le-7 

1  surf  air  shx  shy  shz 
1  ave/surf  1  100000  100000  cl 


57 


dump 

1  surf  100000  dump.vane50rvbb  id  vlx  vly  viz  v2x  v2y 

v2z  v3x  v3y  v3z  f_l 

compute 

2  surf  air  shx  shy  shz 

fix 

2  ave/surf  1  10000  10000  c_2 

dump 

2  surf  10000  dump.vane50tenthrvbb  id  vlx  vly  viz  v2x 

v2y  v2z  v3x  v3y  v3z  f_2 

stats 

1000 

stats_style 

step 

run 

100000 
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